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Abstract 

The QCD analysis of deep inelastic diffractive scattering at HERA is performed assuming the 
dominance of the ’’soft” pomeron exchange and simple, physically motivated parametrization of 
parton distributions in the pomeron. Both the LO and NLO approximations are considered and 
the theoretical predictions concerning the quantity R = Fj^ /for diffractive structure functions 
are presented. 

Resume 

Une analyse de la diffusion diffractive profondement inelastique a HERA est presentee, basee sur 
la dominance d'un pomeron ”doux” et une simple, physiquement motivee parametrisation des 
distributions des partons dans le pomeron. Les approximations des logarithmes dominants et au- 
dela des logarithmes dominants ont etees considerees et les predictions theoriques pour la quantite 
R = Fj^/F^ pour les fonctions diffractives de structure sont presentes. 


The main aim of this paper is to analyse the 
diffractive processes in deep inelastic ep scattering at 
HERA 1^, H, 1^, assuming that they are dominated by 
a ’’soft” pomeron exchange with the pomeron being 
described as a Regge pole with its trajectory 

ap{t) = ap{0) + a'p t , (1) 

where ap(0) w 1.08 and a'p = 0.25 GeV“^0-[|^ . 
Different description of these processes based on the 
’’hard” pomeron which follows from perturbative QCD 
has been discussed in |l^-Q. 

The diffractive structure functions have the follow¬ 
ing factorizable expression : 


dF^L^xp,l3,Q‘^,t) 

dxpdt 


= fixp,t) Flp{(3,Q‘^,t) . 


( 2 ) 


The function f{xp,t) is the so called ’’pomeron flux 
factor” which, if the diffractively recoiled system is a 
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single proton, has the following form [l^ : 

(3) 

Idtt 

where B{t) describes the pomeron coupling to the 
proton and is parametrized as below 0 

B{t) = 4.6 exp (1.9GeV“^ t) . (4) 

The normalization factor N was set to be equal to 
^ following the convention of refs. j|, In 

addition we have slightly increased ap(0) in the flux 
factor f{xp,t) (^) up to ap(0) = 1.1 Q. 

The functions F^piPjQ^jt) are the pomeron 
structure functions with the variable /3 playing the role 
of the Bjorken scaling variable for the 7* *(Q^) pomeron 
inelastic ’’scattering”. 

In the region of large the pomeron structure 
function F^ is expected to be described in terms of 
the QCD improved parton model and is related in the 




conventional way to the quark distributions qf (/3, Q'^,t) 
in the pomeron: 

F[{f3,Q^t) = 2 /3 qfif3,Q\t) , (5) 

i 

where are the quark charges (note that qf = qf ■ ) 
The above formula is given in the leading logarithmic 
approximation of perturbative QCD in which the 
pomeron longitudinal structure function = 0. 

At first we shall specify the details of the parton 
distributions in the pomeron at the reference scale Qq = 
4GeV^. 

At small /3 both the quark and gluon distributions 
are assumed to be dominated by the pomeron exchange 


Pqfi/3,Qlt) = af(<) 

P = a^(t)(6) 

The functions af (t) and {t) can be estimated from 
the factorization of pomeron couplings 

of W = rit) a, Gg (t) = r{t) Gg , (7) 


where the parameters and Gg are the pomeron 
couplings controlling the normalization of the small x 
behaviour of the sea quark and gluon distributions in 

the proton i.e. 


xqi{x,Ql) +xqi{x,Ql) 
xg{x,Ql) 

and the function r(t) is: 


2 a, 

Gg ( 8 ) 


r{t) 


TT Gppp(t) 

2 B{Q) 


(9) 


The coupling Gppp{t) is the triple pomeron coupling 
and its magnitude can be estimated from the cross- 
section of the diffractive production p + p p + X in 
the limit of large mass Mx of the diffractively produced 
system X. We neglected the (weak) t dependence of the 
function r{t) and have estimated its magintude from 
the recent Tevatron data as r{t) w r(0) = 0.089. 
The parameters ai were estimated assuming that the sea 
quark distributions in the proton can be parametrized 
as: 


xqi{x, QI) + xqi{x, QI) =2 g^ x^ (1 - x^ ( 10 ) 


and fixing the constants from the requirement that 
the average momentum fraction which corresponds to 
those distributions is the same as that which follows 
from the recent parametrization of parton distributions 
in the proton 00- The momentum sum rule has also 
been used to fix the parameter Gg i.e. we assumed 

xg{x, Qo) = Og (1 - xf (11) 


and imposed the condition that the gluons carry 
1/2 momentum of the proton. We extrapolated the 
pomeron dominated quark and gluon distributions in 
the pomeron (see (^) to the region of arbitrary values 
of P by multiplying the factor by 1 — /3 0. 

We have also included the term proportional to /3(1 — 
P) in both the quark and gluon distributions j^. The 
normalization of this term in the quark distributions 
has been estimated in assuming that it is dominated 
by the quark-box diagram with the non-perturbative 
couplings of pomeron to quarks. In this model one gets: 

Gtt 

Pq^{P.Ql) = — P{l-P) , (12) 

where C ~ 0.17 |^. We found that the fairly reasonable 
description of data can be achieved provided that the 
constant G is enhanced by a factor equal to 1.5. We 
have also assumed that the relative normalization of 
the quark distributions in the pomeron corresponding 
to different flavours is the same as that of the sea 
quark distributions in the proton 0,0- Finally the 
normalization of the term proportional to /3(1 —/3) in the 
gluon distribution in the pomeron has been obtained 
by imposing the momentum sum rule. Following the 
approximations discussed above we have neglected the t 
dependence in those parton distributions. 

As the result of the estimates and extrapolations 
discussed above the parametrization of parton distribu¬ 
tions in the pomeron at the reference scale Qq = 4GeV^ 
looks as follows: 


P 9^iP.Ql) 
PuP{P,Ql) 
PdP{P,Ql) 
Ps^iP.Ql) 
Pc^iP.Ql) 


(0.218 -b 3.30 P) (1-/3) 

0.4 (1-^) S'^(/3) 

0.4 (1 - ^) S^{P) 

0.2 (1-5) S^{P) 

5 S^iP) , (13) 


where the function S^{P) is parametrized as below 


= (0.0528 /3"° °® + 0.801 /3) (1 - /3) (14) 


and 5=0.02 ||^, The analysis of the pomeron 

structure functions based on different parametrizations 
of parton distributions in the pomeron has recently been 
presented in refs. |p|, 0. 

The parton distributions defined by eqs. ( 00 ) were 
next evolved up to the values of for which the data 
exist using the LO Altarelli-Parisi evolution equations 
10 , 0 with A = 0.255 GeV [|^. In Fig.l we show our 
results for the quantity: 


pXpH /*0 

F^{P.Q^)= dxp 

J xpr J —( 


dt 


dFP{xp,P,Q'^,t) 


dxpdt 


(15) 


with with xpL = 0.0003 and xph = 0.05 plotted as 
the function of xp for different values of /3 and . We 
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Figure 1. Theoretical predictions for diffractive F^{l3,Q^) and 
their comparison with the data from HERA (solid lines). 

Dashed lines show predictions when there are no gluons at the 
initial scale Qq. 


see a very good agreement with the recent data from HI 
collaboration at HERA |^. Dashed lines in Fig.l show 
Q^) computed when gluons were neglected at the 
initial scale Qg = 4GeV^. 

We have also performed the NLO QCD analysis in 
order to be able to consistently introduce the quan¬ 
tity R = FE{f},Q^)/F^{P,Q^) where FE{(3,QE = 
FE{(3,QE - FE{P,Q^) and Ef (/ 3 ,Q 2 ) jg defined in 
analogy to (H). We compared the NLO and LO results 
for FE{P,Q'^) from Fig.l and found that both approx¬ 
imations lead to similar predictions. In Fig.2 we show 
our results for R. The longitudinal diffractive structure 
function is driven mainly by the gluon distribution in 
the pomeron and a large amount of gluons (see im¬ 
plies that R can reach 0.5 for f3 < 0.1. 
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Figure 2. Theoretical predictions for diffractive R plotted as the 
function of (for fixed (3) and the function of 0 (for fixed Q'^). 
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